Critical assessment of two-qubit post-Markovian master equations 



(N 

o 






S. Campbell 123 , A. Smirne 4 , L. Mazzola 1 ' 5 , N. Lo Gullo 2 , B. Vacchini 4 , Th. Busch 2 ' 3 , and M. Paternostro 1 

1 Centre for Theoretical Atomic, Molecular, and Optical Physics, 

School of Mathematics and Physics, Queen's University, Belfast BT7 INN, United Kingdom 

'Department of Physics, University College Cork, Republic of Ireland 

3 Quantum Systems Unit, Okinawa Institute of Science and Technology, Okinawa, Japan 

4 Dipartimento di Fisica, Universitd degli Studi di Milano, Via Celoria 16, 

1-20133 Milan, Italy & INFN, Sezione di Milano, Via Celoria 16, 1-20133 Milan, Italy 

5 Turku Centre for Quantum Physics, Department of Physics and Astronomy, 

University of Turku, FI-20014 Turun yliopisto, Finland 

(Dated: March 1, 2013) 

A post-Markovian master equation has been recently proposed as a tool to describe the evolution of a system 
coupled to a memory-keeping environment [A. Shabani and D. A. Lidar, Phys. Rev. A 71, 020101 (R) (2005)]. 
For a single qubit affected by appropriately chosen environmental conditions, the corresponding dynamics is 
always legitimate and physical. Here we extend such situation to the case of two qubits, only one of which ex- 
periences the environmental effects. We show how, despite the innocence of such an extension, the introduction 
of the second qubit should be done cum grano salis to avoid consequences such as the breaking of the positivity 
of the associated dynamical map. This hints at the necessity of using care when adopting phenomenologically 
derived models for evolutions occurring outside the Markovian framework. 

PACS numbers: 03.65.Yz, 03.67.-a, 42.50.Lc 
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I. INTRODUCTION 

The study of open quantum systems, which has been ex- 
tensive over the years yj y|, is interesting for many stimu- 
lating reasons. First, the occurrence of uncontrolled system- 
environment interactions constitutes a fundamental obstacle 
for the realization of reliable control over quantum devices, 
and the challenge is to design practical schemes for quench- 
ing such unwanted couplings. Second, more fundamentally, 
the dynamics of open quantum systems offer a way to access 
and explore the way genuine quantum features are smeared 
out into classical ones B- 

Very often, the formal analysis of an open quantum dy- 
namics is performed by invoking the use of two simplifica- 
tions: the weak system-environment coupling and the forget- 
ful nature of the environmental system |@, a- This defines 
the so-called Markovian framework, which is often useful 
for the grasping of a qualitative understanding of a system- 
environment evolution and sometimes even physically justi- 
fied by particularly favorable working conditions. Notwith- 
standing its pragmatic handiness, it should be kept in mind 
that Markovianity is only an approximation and the condi- 
tions for its application frequently do not match the reality 
of a given physical situation |6]. The increasing awareness 
of the limitations in the Markovian framework, the identifica- 
tion of explicit cases of non-Markovian system-environment 
evolution and the ability to experimentally simulate structured 
reservoirs with inherent non-trivial dynamics [7] have trig- 
gered the study of open quantum systems beyond the Marko- 
vian regime [8]. A considerable number of techniques, both 
analytic and numerical, have been put forward with the goal of 
accurately tackle non-Markovian evolution II U I9l4l2ll . Among 
them, post-Markovian (postM) master equations have been 
proposed as a means to interpolate between fully Markovian 
Lindblad-like approaches and the exact Kraus-operator pic- 



ture of the reduced dynamics of a system coupled to its en- 
vironment J131U4I1 . offering clear conditions for physicality 
of the dynamics. This circumstance is particularly relevant if 
one considers that the prerequisites for meaningful evolution 
(i.e. the dynamical map evolving the density operator has to 
be trace-preserving, positive and completely positive) are of- 
ten violated by dynamics arising from non-Markovian master 
equation. 

In this paper we show how even a well-tested tool as the 
post-Markovian master equation retains strong limit of appli- 
cability hidden in it. We consider a qubit system embedded 
in an environment whose dynamics is meaningfully described 
by a post-Markovian master equation and an ancillary qubit, 
completely decoupled from the system and the environment, 
whose unitary evolution is ruled by a von Neumann equation. 
As the evolutions of system and ancilla are independent, the 
master equation of the enlarged system containing the system 
and the ancilla can be naturally written as the sum of the gen- 
erators of the two disjointed dynamical evolutions. We show 
how the solution of such a master equation breaks positiv- 
ity and creates correlations between the two non-interacting 
qubits even when initially prepared in a factorized state. Fol- 
lowing the formulation of Shabani and Lidar, we go beyond 
the intuitive approach and derive a post-Markovian master 
equation which accounts ab initio of the ancillary degree of 
freedom. Such a master equation for the extended two-qubit 
systems fails to describe a completely positive or even only 
positive dynamics, as in the previous case. The loss of pos- 
itivity of the two extended evolutions occurs under the same 
conditions for which physical dynamics of the single qubit- 
system is guaranteed [15]. As such our results shed light on 
the limit of validity of post-Markovian master equations and 
warn against the hazard of naive extensions of non-Markovian 
equations of motion. 

The remainder of the paper is organized as follows. Section 
II outlines the postM master equation and its derivation that 



will be focus of this paper. Section III deals with an intuitive 
approach to introducing an ancillary system into the evolution 
and examines the loss of physicality in the map. In Section IV 
we employ the methods of Shibani-Lidar and derive the asso- 
ciated postM master equation explicitly including the ancilla 
in the derivation and show this also fails to give a physical 
map. Section V we discuss the steps necessary to resolve the 
loss of physicality when introducing an ancillary system. Sec- 
tion VI we draw some final conclusions. Finally, there is an 
Appendix in which some of the more technical steps involved 
in the calculations are presented. 



II. POST-MARKOVIAN MASTER EQUATION 

It is instrumental for the second part of our analysis and 
pedagogically quite useful to briefly remind of the key steps 
needed in order to derive the postM master equation proposed 
in H 1411 . The basic observation is that the operator-sum rep- 
resentation of a system-environment dynamics corresponds 
to making a single measurement on the environment at time 
t. In such a measurement picture, a master equation in the 
Lindblad form can be seen as a continuous series of mea- 
surements performed at infinitely closely spaced times [14]. 
The postM approach interpolates such extremal cases: At 
time t = 0, system and environment start their joint evolu- 
tion until, at a random time t'<t, a measurement on the bath 
is performed and the state of the system becomes A(f')p(0), 
where A(f') is a one-parameter map independent of the final 
time t. The measurement resets the state of the environment 
while the joint system-environment evolution proceeds from 
then until time t, when the final measurement is performed. 
The probability for the extra environmental measurement at 
t' to occur depends on the memory properties of the envi- 
ronment itself, which are accounted for by the introduction 
of a memory kernel k(t',t). The final state of the system is 
thus obtained averaging out over the different random times 

as P (t) = ^random t> W , t)K(t - t')p(t'). 

The derivation proceeds with the discretization of the inter- 
val of time [0, t] into segments of equal length e. By evaluat- 
ing the difference of the density matrix operator at two con- 
secutive instant of times, dividing for e and performing the 
limit for e— >0 we obtain the master equation 

p= J dt' [k(t', t)A(t')d (l - n + d,k(t', t)A(t')] p(t - f). 
Jo 

In the spirit of perturbation theory, the relation p(t - t') — 
A(f - t')p(Q) is used inside the integral (meaning that in the 
first order the map describing the dynamics belongs to a one- 
parameter family). Assuming that the inverse of the map 
A(f - t') exists, one can write p(0) = A _1 (f - t')p(t - t') and, 
with the assumption that the memory kernel function depends 
on one parameter only k(t',t) - k(t'), the master equation 
above becomes 



In order to obtain a master equation for each particular prob- 
lem, one needs to specify the form of the map A(f). In the case 
of a Markovian super-operator A(f') = e^' the postM master 
equation attains the form 



dp(t) 
dt 



= t\ dt'k(t')e £t ' p(t-t') 
Jo 



(2) 



dp(t) 
dt 



J dt' 

Jo 



k(t')Mt') 



dMf- 



d(t - f) 
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opit-t'). (i) 



We consider here a qubit interacting in a dissipative fashion 
with a non-zero temperature bosonic bath, the Markovian Li- 
ouvillian £. being in the form 

£p=yn[& + pfr---{fr-&+, p}] +y(n + 1)[<t_p<t + --{<x + <t_, p}] 

(3) 
where & ± are the spin raising and lowering operators. The 
memory function of the bath k(t) is taken of exponential form 
k(t) — X ex P(~X { )- The map O(f) that arises from Eqs. (fj) 
with such a choice of the Liouvillian super-operator was con- 
sidered extensively in Refs. [15] and was demonstrated to be 
completely positive for all range of parameters. Therefore the 
master equation in Eq. ((2j gives a physically meaningful evo- 
lution for a single damped qubit. 



III. INTRODUCING AN ANCILLARY SYSTEM : 
INTUITIVE APPROACH 

Here we consider the extended system comprising a 
damped qubit, whose evolution follows the postM master 
equation in Eq. (O, and an ancillary qubit evolving in time 
unitarily and independently from the first qubit. The intro- 
duction of ancillary degrees of freedom is useful in many ap- 
plications, for example in view of quantifying the degree of 
non-Markovianity of the dynamics of the system [16] (this 
can be done by monitoring the evolution of the entanglement 
between system and ancilla) and for inferring the specifics of 
the system-environment interaction [17] (one can access the 
property of an environment affecting an inaccessible system 
by means of continuous measurements on an ancilla). 

To derive the master equation describing the extended sys- 
tem we use the following line of reasoning. The starting 
point of a microscopic approach is the von Neumann equa- 
tion for the density operator of system, ancilla and environ- 
ment p SAE = -i [h S e ®h A +H A ® i SE , Psae\ containing a 
Hamiltonian term coupling the system and the environment 
Hse, from which decoherence of the system arises, and an in- 
dependent ancillary Hamiltonian term Ha. By tracing out the 
environmental degrees of freedom one obtains 

p SA = -/Tr £ \H S E ® 1a, Psae] - i [h a ® is . Psa] , 

meaning that the generator of the environment-induced sys- 
tem evolution and the one of the ancillary evolution simply 
sum up. As the postM generator accounts for the description 
of the dynamics of the single damped qubit, the most natural 
way to construct the master equation for the two qubits is to 
substitute the term corresponding to the dissipative dynamics 




FIG. 1: (Colors online) (a) Trace distance T (dashed line) and lowest eigenvalue E/ (solid line) of the mapped matrix pit) solution of Eq. <(4j 
against the rescaled interaction time t = Jt for the initial state p(0) = \0} s ® (|0> + |1)W v2 and 7 = 2. (b) Lowest eigenvalue e; of the matrix 
pit) as a function of rescaled time t = yt and the ratio K = J/y for the same initial state of panel (a), (c) Plot of the lowest eigenvalue £ rf of 
the Choi matrix associated with the dynamics of pit), (d) Lowest eigenvalue si of the matrix pit) solving Eq. 0. In all the panels we have 
taken n = 1 with;f=y. 



with the postM generator of Eq. (|2}, thus obtaining 
dp(t) 



dt 



-=-l[i/,p(0]+£s®lA 



, fdt'k(t')e ls ''<i 
Jo 



UfKt-f), (4) 



with H = Ha® Is and pit) the density matrix for the system- 
ancilla pair. A first check supporting the reasonableness of 
the previous equation is that when tracing out over the ancilla 
(system) degree of freedom we obtain the postM equation of 
the damped system qubit (von Neumann equation for the an- 
cilla). In order to be meaningful this master equation has to be 
trace preserving and has to give rise to a physical dynamics, 
meaning that p(t) has to be a density matrix (trace one, Her- 
mitian and positive) at any time and the map associated has to 
be completely positive. We choose the following form for the 
ancillary Hamiltonian 



H A = J&f, 



(5) 



where J is the frequency split between the logical states of the 
ancilla and <x^ is the corresponding z-Pauli operator. In order 
to solve Eq. (0|l we decompose p(t) in terms of the damping 
basis 1 18] of X, as 



Pit) = Yj C 'i< f )0-f 



(W = o,z,±) 



(6) 



',] 



where &.' are the right eigen-operators of £$ such that 

£ s &f = Ai&f [with o-o = (l/2)(i - ctj{2n + 1))] and c u (t) 
are time-dependent coefficients. The choice of such an op- 
erator basis for the density matrix is very convenient, as the 
action of X.s on each element of the decomposition of pit) 
simply depends on the eigenvalues associated with the damp- 
ing operators. Moreover, it turns out that the Hamiltonian of 
Eq. (O is diagonal with respect to the damping basis. By in- 
troducing the dual of the damping basis [18] (defined as the 
solution of the orthonormality condition Tr[<r,<r ; ] = 6tj, with 
6/j the Kronecker delta function) it is possible to obtain a set of 
16 equations for the time-dependent coefficients c,j{t). These 
equations can then be solved by using Laplace transform tech- 
niques. Once we have at hand the solution of the master equa- 
tion, we can analyze its physicality, i.e. whether the associ- 
ated map is completely positive and p(t) is a density operator 



at all times. Moreover, it is interesting to compare the solu- 
tion of this master equation p(t) for initially factorized states 
with the tensor product of the two independent system-ancilla 
evolutions psif) ® PAif), obtained solving the postM master 
equation for the system in Eq. (O and von Neumann equation 
for the ancilla pA(i)= —i H A , Pa(0 L respectively. To this aim 
we make use of the trace distance Tip, (p)=^Ti[\p - (p\] (here 

\A\- VA^A is the trace norm of an arbitrary matrix A), which 
defines a distance between the density matrices p and ip. 

We start our analysis by considering the system-ancilla ini- 
tial state p(0) = |0> s ® (|0> + |1>W VI Figure[T](a) shows 
the trace distance (dashed line) between pit), the solution of 
Eq. (H|, and p s (t) ® Pa(0 an d tne ti me evolution of the small- 
est eigenvalue of pit) (solid line), for the set of parameters 
n=l, x-J-JI^-- As we can see, not only pit) differs substan- 
tially from the meaningful expected solution but also violates 
positivity, being its smallest eigenvalue always negative. This 
rather surprising result implies that the matrix pit) ceases to 
be a density matrix and that Eq. (0]i fails to describe a physical 
dynamics. As the complete positivity of the map of the sin- 
gle damped qubit (derived by Eq. (01) is guaranteed for any 
values of the parameters n, y and %, we ascribe such a dra- 
matic implication to the introduction of the ancillary term and 
check whether there is a regime of the J parameter in which 
physicality is saved. To this purpose we investigate the evo- 
lution of the same system-ancilla factorized initial state as a 
function of the J parameter. Figure [2(b), showing the lowest 
eigenvalue of the matrix pit), demonstrates that positivity is 
broken even for small value of J. This result leaves no doubt 
about the lack of physicality of the evolution described by the 
Eq. ©. For the sake of completeness we include the check of 
the (more restrictive) complete positivity condition. In order 
to test if a map is CP one can use Choi's theorem [ 19] which 
states that a map <& is CP if the matrix Co = Z,y Ejj®<$iEjj) is 
positive, where Ey is a matrix with 1 in the //'' entry and ze- 
ros everywhere else. The map O can be determined from the 
analytical solution (see Appendix). Figure [T](c) presents the 
lowest eigenvalue arising from Co as a function of time and J, 
and confirms that the map loses complete positivity whenever 
J + 0, as expected. 

We notice that the evolution map associated to Eq. (0|i is not 



given by the tensor product of two individual evolution maps 
of the system and the ancilla, as one would have hoped for. 
Indeed the overall evolution map is not obtained as the com- 
position of the completely positive maps individually gener- 
ated by the two terms at the right hand side of Eq. ® , despite 
the fact that they do commute and have a common set of time- 
independent eigenoperators. In the attempt to trace back to 
the reasons of the break of positivity and gain insight on such 
a fictitious evolution, we can calculate the local in time gen- 
erator of the dynamics. This corresponds to consider the op- 
erator pit) as the solution of the equation pit) - K TCL it)pit), 
therefore since we have that pit) - <D(f)p(0), the local in time 
generator will be equal to KjcL(t) - <J>(t)Q>(t)~ l . The expres- 
sion of KjcLit), presented in the Appendix demonstrates how 
the addition of the commutator with the ancillary Hamilto- 
nian leads to the appearance of terms (of the form df ® <r*) 
correlating the two independent qubits. In this light the trace 
distance plotted in Fig. [T](a) not on ly describes how far the 
evolution map O(f) is from its factorized counterpart, but also 
accounts for such unphysical correlations created between the 
qubits. 



after tracing over the system degree of freedom, one obtains 
p A = - i [H A ,^dt'kit')e- if,A, 'p A it - t')e iAA '']. This equation 
cannot be cast in the form of a von Neumann equation, there- 
fore it cannot describe the unitary evolution of the ancillary 
qubit. Furthermore, such equation does not even describe a 
physical evolution, since it breaks positivity. The reason for 
that can be tracked down in the founding idea of the postM 
approach: the use of an additional measurement on the bath 
at a random time perturbing the dynamics described by the 
exact Kraus operator representation. Even if the ancillary sys- 
tem does not interact with the bath, the randomicity of such 
an instant of time (weighted for the memory function of the 
bath) modifies the purely unitary dynamics of the ancilla, no 
longer recovered when tracing over the system. This analysis 
demonstrates that the verbatim application of the postM recipe 
to systems where a partially unitary dynamics is involved risks 
to produce highly unphysical results. 



V. INTRODUCING AN ANCILLARY SYSTEM: CURING 
PHYSICALITY 



IV. INTRODUCING AN ANCILLARY SYSTEM: 
SHABANI-LIDAR APPROACH 

As the system dissipator and the ancillary unitary term act 
on two different Hilbert spaces, the addition of the two inde- 
pendent generators appears as the most natural way to build 
the master equation of the extended system. Nevertheless, 
one may ask whether the results obtained are consequences 
of a possibly naive approach, and an alternative description 
can cure the physicality of the evolution. Here we include 
the ancillary system and its Hamiltonian evolution among the 
building blocks of the construction of the postM master equa- 
tion. 

We take Eq. (fTJ and consider A as the composition of the 
Markovian dissipative generator As it) = e Ls> acting on the 
system and the unitary map K A it)ip) = U A (t)(p)Ul(f), with 
U A it) - e~' HA ' generated by the Hamiltonian H A . With this 
substitution and performing the necessary calculations, we ob- 
tain the following master equation 



—pit) = £ s ®h 
at 



a Zst'-iH A t' 



)e J - s ' e-'" A ' pit-t')e 



i\JH A t' 



-I 



I dt'kif 
Jo 

i s ®H A , \ dt'kit')e tst 'e~ mk< pit- t')e iHAf 
Jo 



(7) 



We start checking the positivity of the evolution of the initial 
state p(0) = |0) s ®(|0) + |1»a/ V2 as a function of the J param- 
eter. The plot in Fig.Q~|(d) shows that positivity is broken also 
in this case. However, there is an even more important pre- 
liminary physicality check that fails to be satisfied: the trace 
over the system degrees of freedom of Eq. (IT) does not pro- 
duce the von Neumann equation for the ancilla. To illustrate 
this statement we consider the simple case where y — in the 
dissipator X.s ■ For this choice of parameter the first term at 
right hand side of the master equation in Eq. (0 is null and, 



We finally show why the naive generalization of the postM 
master equation failed. Let us consider Eq. ([TJ where now 
pit) is the density matrix of the system-ancilla pair in the in- 
teraction picture with respect the bath (Hb) and system-ancilla 
(Ho) free Hamiltonians. We take A(f) = A s (t) ® t A with, as 
before, As(f) = e^ s '. By performing a transformation to the 
Schrodinger picture by means of the operator f Uo(f) - e~' H °' 
and assuming that system and ancilla are not interacting and 
their initial state is factorized, we get 



pit) = -i[H s ,p s it)] ® p A {t) - ipsit) ® [H A , p A it)] 

+ ls [ dt'k(t')e lst, €l s (t')p s (t - f)€ll(t') ® p A it), 
Jo 



(8) 



where, with a little abuse of notation, the density matrix pit) is 
meant to be in the Shrodinger picture, Hsj A ) is the free Hamil- 
tonian of the system (ancilla) and f Us(t) - e~' Hs ' . The pres- 
ence of the term e^ s '' f U s (t')psit - f)€(l(t') tells us that after 
the extra generalized measurement preformed by the bath at 
time t' the evolution of the system proceeds for a time interval 
t—f under the Markovian map [ 14]. The above equation thus 
generates a global map O(f) = O s (f) ® A (f) where O s (f) is 
the one induced by the master equation (f2]i (in the Schrodinger 
picture) and <b A (t)p A -> €(oit)p A €(lit). The positivity of O(0 
is guaranteed by the positivity of the single maps once con- 
ditions for complete positivity of Os(f) are met 11411 . <J>(f) is 
thus a proper generalization of the postM master equation as 
one would expect. 

In order to derive Eq. (|8) ab initio, one has to propagate the 
effects of the bath from time f up to t UM . This obviously 
makes sense only for those degrees of freedom which are ac- 
tually interacting with the bath. This is not the case for the 
ancilla, whose dynamics is entirely unitary. It is then not suf- 
ficient that the super operator A(f) does only act over the rele- 
vant part of the Hilbert space (A(f) = Aj(f) ® l A ): the propa- 



gation has to involve only that particular subspace. This con- 
dition cannot be met if one derives Eq. (O in the Schrodinger 
picture by simply replacing A(f) , given that the memory ker- 
nel "keeps track" of all of the degrees of freedom regardless 
wether or not they did interact with the bath at one point in the 
past. This problem can be overcome by introducing a memory 
kernel which is of the form k{t)tM ® §(f)tp where 1m is the 
identity operator over the subspace the bath "has memory of" 
whereas 1^ does act over the "free" subspace. 



VI. CONCLUSION AND FINAL REMARKS 

In this paper we highlighted the hazard of using postM mas- 
ter equations in systems containing a partially unitary dynam- 
ics. The specific example of a physical postM master equa- 
tion describing a damped qubit extended to include an ancil- 
lary one shows that the introduction of a Hamiltonian term 
has drastic effects even when the ancilla is unitarily evolv- 
ing and totally decoupled from the system and the environ- 
ment. Both the intuitive approach, based on the summation of 
the two independent generators of system and ancilla, and the 
postM derivation, accounting for the ancilla degree of freedom 
ab initio, break fundamental physicality requirements such as 
positivity. As final remarks we add that the same conceptual 
problems arise when using Hamiltonian terms coupling sys- 
tem and ancilla (such as Ising coupling), and also in the case 
where no ancilla is brought in and the Hamiltonian term acts 
on the system only. Analogous results can be found when the 
equation to generalize is a phenomenological memory-kernel 
master equation. This difficulty has also been considered in 
a different context in [20], where complete positivity for an 
initially well defined integrodifferential master equation has 
been preserved upon the introduction of a Hamiltonian term 
only by suitably modifying the integral kernel. Due to their 
very construction, those master equations do not always guar- 
antee positivity and CP, even in the single-qubit formulation. 

We have addressed the use of postM master equations in the 
case of composite (interacting and non-interacting) systems, 
pointing out that quite a careful approach should be taken to 
include the unitary part of the evolution in the measurement- 
based picture at the basis of the postM framework. For the 
case of non-interacting particles, only one of which experienc- 
ing the effects of the postM environment, the intuitive combi- 
nation of the separate dynamics of the two systems may lead 
up to unphysical results. Our formal finding, which are inter- 
esting from a fundamental viewpoint, will have implications 
at the pragmatic side in schemes where the ancilla is used in 
order to probe the properties of the the environment into whi ch 
the system is immersed [17]. 
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APPENDIX A - Derivation of the time-local generator 

Here we present the derivation of the time-local genera- 
tor associated to the master equation in Eq. (4). By us- 
ing the decomposition of p(f) in terms of the damping ba- 
sis {<x s ) (/ = 0, 3,± with (x^ the z Pauli matrix) given by 
p(t) = ^ijCij(f)&? ® & A > we can derive the following set of 
differential equations 

c u (t) = -i J] ft»GH(0 + k f dTk(T)e A - T c u (t - t) (A-l) 



and <t, ' the elements 



with hfj = Tr {<rf ® & A [h, df ® frf ] 

of the dual damping basis for the system and ancilla. These 
equations can be solved in the Lapace space as 



Cij(t) = Lap 



s + ihj--^ 



Cij(0). 



(A-2) 



where, according to H as given in Eq. (5), one has 
hf; -'■ 5ik6jihj and the time-dependence of the Laplace anti- 
transform Lap 1 [■] is understood. The solution shows that the 
evolution map is diagonal in the basis l&f <S> <r A ) and can be 
written as a diagonal 16x16 matrix 



<D(f) = diag{£ 7 (f)} 



with elements 



&(0 = e 



-lty/2 



cosh - 



/A;;? 



\ 2 



a~ 



■y • , i A ij l 
+ — - sinh ■ 



(A-3) 



(A-4) 



where A, 7 = J(x + At) 2 - h 2 j - lihfx - Ad and £l ij =x-A i +ih j . 

Evidently, the coefficients £y(f) cannot be written as £i(0£/(0> 
therefore the evolution map cD(f) is not given by the tensor 
product Qk^sA^kit) °f single-qubit maps 0^(0- Having the 
matrix representation of the evolution map, we can calculate 
the time-local generator ^tcl(0 = <$>(f)G>(t)~ , whose matrix 
representation in the two-qubit damping basis is given by 



K(f) = <b(f)<J>(0 = diag 



&(0 

fy(0 



(A-5) 



Given a basis {B fl ! , in £.(f{ N ), one can define the 

scalar product between two super-operators defined in the 
space of linear operators 2{"H N ) ||5|] as 



<A!,A 2 > = £ Tr (A! [£„]%[£„]}. (A . 



6) 



A basis in the space of super-operators is given by the family 
[fay} B =0 N 2_ x such ft 31 VjSylp] = BpgB], with g a generic 



density matrix. Therefore, any super-operator A acting on 
&CH N ) can be written as 



fir 



Py 



with Ap 7 ={ipp 7 ,A). We take a basis JB fl } 2 such that 

B = -4=1 and Tr{fi fl ) = 0, Va = 1 , . . . , AT 2 - 1. From the 
expression 



KjcdOlp] = J] Kp y {t)B p pBl (A-8) 

,8,7=0 

of the time-local generator, we obtain its Lindblad-like ver- 
sion simply by removing the first row and the first column of 
the corresponding matrix of coefficients, since the evolution 
map is trace -preserving and thus Tr{KTci,(t)[p]} - 0. We thus 
obtain 



N 2 -l , 

K TCh (t)[p] = -i [H',p] + J] Kpy(i)\Bp P B\ - - [B\Bp,p) 



with the Hamiltonian term given by H' - i-.(g 1 " - g) and 



(A-9) 



N 2 -l 



=— y 



Kpo(f)Bp. 



(A-10) 



In the basis given by 



fig = 1/2 

B s a A = (i S ®df)/2 ; = 1,2,3; o= 1,2,3 






(A-ll) 



7^ A = (df ® 1 )/2 j=l,2,3;a = 4,5,6 

fl^ = (&f ®&fi/2 i,j = 1,2,3; a = 7, ...,15, 



where <Tj ^ are the x and y Pauli matrices respectively, the 
local generator in Lindblad form is given by Eq. (1A-91 with 

K firW = 2 Tr (^l^TCLtAj} (A-12) 



with ^TCL(0[^a] that can be explicitly calculated using the 
relation 

KtclIpI = V y^rTr [&f ® & A p) &f ® & A (i, y = 0, z, ±). 

(A- 13) 
The Hamiltonian term is obtained substituting the expres- 
sion Kpo(t) = i y,„ Tr (BlB]g T all4 in Eq. dA40t. Due 



to the particular choice of the basis in the space of the super- 
operators, the dynamics can create correlations between the 
two qubits only if there are non-zero coefficients Kp y for 
B, y - 7, . . . , 15, i.e. if the dynamics involves basis elements 
of the form &f ® cr A . Thus, a necessary condition for the gen- 
eration of correlations is that there are non zero values outside 
the 7x7 top-left corner reduction of the coefficients matrix 
given by Kp y . For example, one can directly check that if 
^j(t) = £i(t)£j(f) then K Py = for B,y = 7, . . . 15, while this 
is not the case for £/?) given by Eq. dA-4b . In this case, there 
are non-zero coefficients related to & s ® cr A in the dissipative 
part of the generator and there is also a non zero coefficient 
associated to &? ® & A in the Hamiltonian part. This means 
that the Hamiltonian term added in the master equation given 
by Eq. (4) modifies entirely the structure of the time-local 
generator. 



APPENDIX B - Ab initio extension of the post-Markovian 
master equation 

Here we present the mathematical calculations we carry 
over to derive the Master Equation in Eq. (7), obtained by 
including ab initio the ancillary system and its Hamiltonian 
evolution. We recall the notation used: the Markovian dissi- 
pative generator is As(t) = e^', while the unitary map acting 
on the ancilla is A A (t)[p] - UA(f)[p]U\(t). The starting point 
is Eq. (1), namely 



— = \ dt'k(t' 
dt Jo 



)A(f) 



dA(t-t') 



d(t - f) 



-^-K-\t-t')p{t-t'). (B-l) 



As the system and ancilla are not interacting we set A(t) - 
A s (t) ® A A (t). We can thus write 

A-\t-t')p(t-t') = A s \t-t')Ul(t-t')p(t-t')U A (t-t'), (B-2) 

and thus, given that ^g=p = £ s e £s<f - , ' ) = -Cs A s (t - f), we 
have that 

A{t-t')A- x {t-f)p{t-t') = £ s p{t-t')-i[H A ,p{t-f)]. (B-3) 

Using the previous results we finally obtain 

d f 

—pit) = £ s ®h dt'k(t')e £s ' e-' HA ' p(t - t')e lHA < 
dt Jo 

-i i s ®H A , \ dt'k(t')e ls, 'e- iHA, 'p(t-t')e iHA '' . 
Jo 

(B-4) 
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